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Abstract 

A quantum computer is a multi-particle interferometer that comprises beam split- 
ters at both ends and arms, where the n two-level particles undergo the interactions 
among them. The arms are designed so that relevant functions required to produce 
a computational result is stored in the phase shifts of the 2 n arms. They can be 
detected by interferometry that allows us to utilize quantum parallelism. Quantum 
algorithms are accountable for what interferometers to be constructed to compute 
particular problems. A standard formalism for constructing the arms has been devel- 
oped by the extension of classical reversible gate arrays. By its nature of sequential 
applications of logic operations, the required number of gates increases exponentially 
as the problem size grows. This may cause a crucial obstacle to perform a quantum 
computation within a limited decoherence time. We propose a direct and concurrent 
construction of the interferometer arms by one-time evolution of a physical system 
with arbitrary multi-particle interactions. It is inherently quantum mechanical and 
has no classical analogue. Encoding the functions used in Shor's algorithm for prime 
factoring, Grover's algorithm and Deutsch-Jozsa algorithm requires only one-time 
evolution of such a system regardless of the problem size n as opposed to its stan- 
dard sequential counterpart that takes 0(n 3 ), 0(n) and 0(n2 n ). 



A computation entails encoding of a function whether classically or quantum mechan- 
ically. The encoding of a function has been carried out in the form of a bit-flip. Such an 
"oracle" U c in reversible classical computers and also in proposed quantum computers is a 
transformation of an n-bit input x and a work bit w, 

U c : \x) \w) -> |x) \w © f(x)) , (1) 

where © stands for exclusive-OR. When the work bit is initially set to be 0, the oracle 
returns the function value f(x) in the work bit. The customary construction of such an 
oracle pertains to sequential application of reversible gates'. In quantum computation, 
however, the oracle is often converted into a transformation of the form, 

U q : \x) - (-1)'<*> \x) , (2) 

where the information about f(x) is encoded in the phase of a linear superposition state so 
that quantum parallelism could be utilized!@&0. The conversion of U c into U q is performed 
by supplementary transformation! or automatically by appropriately initialized work bit!. 
The construction of the oracle ([I]) by sequential application of one-bit and two-bit gates 
demand an exponential or polynomial number of operations as the problem size grows as 
shown in Table |l|. In order to perform a quantum computation, we need to complete the 
construction of the oracle while the coherence of the system is maintained, although classical 
computation does not bring the issue of a limited decoherence time into a question. We 
propose a concurrent construction of the transformation U q by only one-time evolution of a 
physical system that has arbitrary multi-particle interactions. The exponentially hard work 
that is expressed by the complexity of gate arrays to be applied to a physical system in the 
case of the sequential construction will be replaced by adjustment of an exponentially large 
number of coupling strengths in the system before the coherence of the system is required. 
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The system required to concurrently implement an arbitrary n-bit Boolean function f(x) 
for an n-bit string x consists of n two-level particles with arbitrary multi-particle interactions 
among them. The Hamiltonian of the system is, 

n 

H = Y1 h^iz + H ^ijViz ® a jz + ^ijkViz <8> CT jz (g) O kz 

i=l i<j i<j<k 

H h TuOl2-nVlZ ® °2z® ■ ■ ■ ® CTnz, (3) 



where <S> stands for tensor product. The eigenstates of the Pauli matrix <7j 



(l \ 

for the i-th particle is used as the computational basis {|0) , |1)} for the i-th bit Xj, where 
Ciz = (—l) Xi \xi) (xi = {0,1}). The computational basis of x is defined in the 2 n - 
dimensional state space spanned by the two-basis states of the n particles, \x) = \x n ■ ■ ■ X\) = 
\x n ) ® ■ ■ ■ © \xi). The state \x) represents a number x = Yh=i 2 l_1 Xj. The Hamiltonian is 
diagonal in the computational basis {|00 • • • 0) , |00 • • • 1) , • • • , |11 • • • 1)}. Diagonal elements 
of the (2™ — 1) terms in the Hamiltonian and a 2 n -dimensional vector (1, 1, ■ • • , 1) (diagonal 
elements of the 2 n x 2 n identity matrix) constitute an orthonormal basis to expand 2 n - 
dimensional vectors. It suggests that the transformation U q in (Q), which is a (2 n x 2 n )- 
diagonal matrix, can be constructed by a global phase shift and one-time evolution of the 
system, U-h(t) = e~% HT for time r as 

U q = e^U H (r). (4) 

The concurrent construction of the transformation U q for a given function comprises a prepa- 
ration of the Hamiltonian (by adjusting the coefficients, hu>i, huiij, ■ ■ •, hu>i2... n ) so that (|J) 
is satisfied and a time-evolution by that Hamiltonian. Only the latter process requires the 
conserved quantum coherence of the system. 

The diagonal elements of (H) are decomposed as 



e i7r/0z) _ e i<t> x TT e -WiT{-l)*i x TT^rH)-^ 
i i<j 

i<j<k 

-iuj 12 ... n T(-iri+ x 2+-x 



x ... xe ~Hr— , ( 5 ) 



which gives the solution, 



0= ^£/0)> 

^ a; 

^ = -^E(-ir/(x), 

Z X 

^r = -^£(-ir^/(z), 

= -^E(- 1 ) X ' 1+X2+ - +X '"/(X), (6) 

and determines the coefficients in the Hamiltonian to be prepared so that only one-time 
evolution for time r by the Hamiltonian will construct U q by itself for a given function f(x). 
Note that when f(x) is a constant function, all coefficients in (|6|) except <fi are zero. 

The formalism for concurrently constructing the transformation U q can be applied to 
encode functions and calculated values that are used in existing quantum algorithms. Ex- 
amples are shown below. 

Deutsch-Jozsa algorithn^^M. The algorithm solves the following problem by quantum paral- 
lelism. Given the oracle for an ra-bit Boolean function / : {0, l} n — > {0, 1}, determine either 
(A) / is a constant function (at or 1) or (B) / is a balanced function (the sequence /(0), 
• • •, /(2 n — 1) contains exactly 2 n ~ 1 zeros and 2™ _1 ones). In the original algorithm, the 
oracle that has a form of U c is used twice together with another unitary operation on the 
work bit S : \w) — > (— l) w \w) between the two applications. The three unitary operations 
applied to the system in series are designed so that the information about f(x) is transferred 
from the work bit to the phase of the control register \x) as in (^). Now we have a method 
to construct the transformation U q concurrently, and the work bit can be removed. 

When f(x) takes only or 1, f(x) in (^) can be replaced by —f(x) since 7xf{x) = 
—7if(x) (mod 27r). Therefore a set of parameters obtained by replacing f(x) by 
(-l) Xl+X2+ '" +Xn (f(x) - 2N X ) (N x is an integer) in (D is also solution to (H). For a bal- 
anced function, we choose N x = X1X2, so that co>i 2 ...„r = — ^ (2 n_1 — Y^x^Ng.) = 0. Thus, 
to encode a balanced n-bit Boolean function, we need the multi-particle interactions up to 
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(n— l)-th order (the n-particle interaction is not required). Other parameters are determined 

as, 

^T.(- l T 1+X2+ '" +Xn U^) -2N X ), 

* X 

:. (7) 

Grover's algorithm^. The algorithm explains how a data can be found out of 2 n random 
data entries. The data search problem is described by a function f(x) that returns 1 for 
a single unknown value of x, say x — r, and for the rest of x. The algorithm uses the 
information about f(x) encoded in phase of the control register as in (Q), which can be 
implemented concurrently, when the coefficients in the Hamiltonian are chosen as <fi = 
UiT = -|r(-l) T % WijT = -^(-1)^, • • •, and u 12 ... n r = -£(-l)TL+*H-~r». 
(iii) Simon's algorithm! determines whether a given function / : {0, 1}™ — > {0, l} m with 
m > n is periodic f(x) = f(x') <-* x' = x © s (with a nontrivial string s) or one-to-one. 
Encoding the function f(x) in the oracle of the form U c in the original algorithm can be 
replaced by 

\ x ) ^ e ia£=rf<*) \ x ) ? ( 8 ) 

which can be constructed concurrently. The transformation (H) has the same form as U q , 
where (—1) = e m is replaced by e^ 2 " 1 1 . The formalism for constructing the transformation 
U q works exactly the same way. In order to implement an arbitrary function in this problem, 
all multi-particle interactions up to n-particle interaction in the Hamiltonian are necessary. 
Shor's algorithm for prime factorizatiorE. In order to factorize an odd number N, we 
randomly choose a (N and a need to be relatively prime) and find the order r of a, the least 
r that satisfies a r = 1 mod N. Finding the order is the prime part of Shor's algorithm and 
starts with the transformation on two n-bit registers, 

^ E |a x mod AO , (9) 



where q = 2 n satisfies iV 2 < q < 2N 2 . Replacing the transformation (§) by, 



5-1 ^ q-l 



y \ x \ _, _!_ y mod iv) i v , 10 x 

which encodes a x mod A/" in the phase of the control register instead of in the work bit, 
leaves the algorithm unchanged^. The phase factors a x mod iV in (|10|) can be calculated 
classically as, 

jm( aX mod N ) = e *sf rC=i(«* 2i_1 mod N T\ (11) 



where products refer to multiplication mod (N). Equation (11) has the same form as (|5|) if 
Af = + ±^(-1)^ (Ai = a 2 ' -1 mod iV) is used. Therefore, the time-evolution of the 
system for time r constructs the transformation (|TT| ) by itself when the coefficients are chosen 
as, = n?=i ^ -^r = x oJijT = x x • • ., and -w 12 ... n r = n™=i ^ 

Quantum Fourier transform. To perform Quantum Fourier transform A g for an n-bit register 
(q = 2 n ), 

A q :\x) ^J2e 2 ™ y/q \y), (12) 
y=o 

as used in Shor's algorithm, we need the Walsh-Hadamard transformation on each bit, 

'i i ^ 

in the computational basis, and controlled-phase-shift operators on pairs 



v 1 - 1 / 



of bits, defined as Sjk = e k ~i XiXk on the j-th bit and fc-th bit with j < k where Qk-j — 
7i /2 k ~ J . The right-hand side of @ is equal to Sj± if we choose 2iOjkT = Ok-j, 2ujjT = —6k-j, 
2ujkT = —Qk-j and the rest of the parameters are zero, aside from the global phase factor. 
A product of multiple controlled-phase-shift operators Sip Sip ^\ ■ ■ ■ 5^+1 is also diagonal in 
the computational basis and can be implemented concurrently by setting 2u>i >m T = Q m ~i 
(m — I + 1, • • • , /'), 2oj m r = —9 m -i, 2ujit = — J2 m ®m-i- One-bit rotations and two-particle 
interactions in the Hamiltonian are sufficient to implement controlled-phase-shift operators 
and products of those concurrently. 

The time evolution of a system that concurrently implements functions by our proposed 
scheme and the number of necessary elementary gates for the function implementation by 
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means of the standard sequential scheme are compared in Table [I] for existing quantum algo- 
rithm. It is challenging to find a system that has multi-particle interactions with reasonably 
large and controllable strengths in order to implement arbitrary functions, but if found, 
many functions and quantum algorithms will be implemented by only one-time evolution of 
the system, which may cross out the current biggest obstacle to quantum computation, the 
short decoherence time of a quantum system. 
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TABLES 





Concurrent 






Sequential 


Algorithm /Function 


implementation 






implementation 












Terms in Evolution 


Number of 




the Hamiltonian 


time 


elementary gates 




U q : \x) 


->( 


_!)/(*) \ x ) 




General Boolean function f(x) 
Deutsch-Jozsa 


RI m . . . /(n) 


0(1) 


0(n2 n ) 


Grover 


RI m . . . /w 






O(n) 




\x) -> exp [^/(x)] | 


x) 


\x)\Q)^\x)\f{x)) 


Shor (ip = 2^r, f(x) = a x mod N) 


mm . . . /w 


0(1 


) 


0(n 3 ) 


Simon (<p = ^£=t) 








0(mn2 n ) 






exp 


TT 


\x) 


controlled-phase-shift Sj^ 


RI& 


0(1) 


0(1) 


H n S n - 


-l,nH n —iS n —2,nS n —2,n—lH n 


-2 • • 


H2Si y7l Si y 


n-l ■ ■ ■ Sl^Si^Hi 


Quantum Fourier Transform 




O(n) 


0(n 2 ) 



Table 1. 
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Scaling laws of evolution time and number of elementary gates necessary to implement 
functions for existing algorithms by means of the proposed concurrent implementation and 
the standard sequential implementation, respectively. In the most general case, an n-bit 
Boolean function f(x) : {0, 1}™ — > {0, 1} is encoded in phases of the control register \x), as 
Uq : \x) — > (— \x), in order to utilize quantum interference effect. To implement such an 
n-bit Boolean function, we prepare the n-particle system that has multi-particle interactions 
among them, from two-particle interactions 1^ up to n-particle interactions 1^ in addition 
to one-bit rotations, R. A one-time evolution by the Hamiltonian with properly chosen 
coefficients of the terms in it builds the transformation U q by itself. A Boolean function 
f(x) used in Deutsch-Jozsa problemS has a constraint that f(x) is either constant (at 
or 1) or balanced. Because of the constraint, the concurrent implementation does not call 
for the n-particle interaction I^ n \ In Grover's data search algorithm!, the function to be 
implemented f(x) is zero for all x but r (/(r) = 1) which we search for. As a general Boolean 
function, this function requires the system that has all multi-particle interactions in order to 
be implemented in a concurrent fashion. The transformation U q is constructed by successive 
applications of elementary gates (one-bit gates and two-bit gates) by the standard sequential 
implementation with the help of a work bit prepared in the state ^(|0) — |1)) and another 
qubitB. The number of required gate operations to implement an n-bit Boolean function in 
the sequential manner scales as 0(n2 n ) as opposed to a one-time evolution of a system in 
the case of the concurrent implementation. The function implemented in Grover's algorithm 
is special in that it only needs one (n + l)-bit gate that is constructed by 0(n) elementary 
gates. In Shor's algorithm for prime factoring an odd number ivi, a x mod N (a is a randomly 
chosen integer relatively prime to N) is encoded in a work bit, as \x) |0) — > \x) \a x mod N). 
The construction of a gate array for this transformation requires 0(n 3 ) elementary gates 
by means of the sequential implementation. Instead of encoding \a x mod N) is the work 
bit, encoding it in the phase of the control register \x), as \x) — > exp i^(a x mod N) \x), 
still wor ks§. The factor ^? is determined to differentiate all possible values \a x ) N takes 
(between and N — 1) and also to avoid a destructive interference that ruins the algorithm. 



In Simon's algorithm!], a function to be implemented is / : {0, 1}" — > {0, l} m (m > n). The 
sequential implementation of such a function in a work bit is in need of 0(mn2 n ) elementary 
gates. As in Shor's algorithm, the function / can be implemented in the phase of the control 
register \x), as \x) — > exp ^hrf( x ) ■ To distinguish all 2 m possible values / takes, the 
phase space of 2tt is divided by 2 m (a Boolean function, / : {0, l} n — > {0, 1}, is a special 
case where m = 1. In both Shor's and Simon's algorithms, the concurrent implementation 
of necessary functions involves all multi-particle interactions I^ 2 \ ■ ■ •, 1^ in the system. A 
controlled-phase-shift operator Sj t k acts on a pair of bits, in this case j-th and k-th qubits 
of the control register \x). It adds a phase factor exp (i-^rj^J only when both Xk and Xj 
are ones. Only two-particle interactions between k-th and j-th particles to implement 
this operator. Such controlled-phase-shift operators on all pairs of qubits (n(n — l)/2 paris 
in total) compose Quantum Fourier transform, which is used in Shor's algorithm, along 
with the Walsh-Hadamard transformation on each qubit (Hi on l-th qubit)!. The Walsh- 
Hadamard transformations are applied to all qubits, from X\ to x n , and controlled-phase-shift 
operators Si jn Si jn -i ■ ■ ■ Sij+x are interleaved between Hi and H + i. A controlled-phase-shift 
operator and a product of multiple of them can be concurrently implemented by a one-time 
evolution of a system that has only two-particle interactions and one-bit rotations. In total, 
Quantum Fourier transform is constructed by n Walsh-Hadamard transformations and n 
time-evolutions of a system in a concurrent manner, whereas in the case of the sequential 
implementation, n(n — l)/2 controlled-phase-shift operators are implemented in series. 
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10. Following Shor's original algorithm!!, We then apply the Quantum Fourier transform 



( [12]) to the control register \x) in the state (fLOD, and then we obtain 

i 9-1 9-1 

_ mod N ) e 27rix y/i \ y ) . 

Q x =0 y=Q 

A measurement projected to \y) has a probability, 

1 9-1 

Since the order of a is r, a x mod N may be written as a k , < k < r, where x = k mod r. 
Therefore, the probability is, 



e i "sw^ a,x mod N ) e 27Tix y/i 



^e'w^ mod N )- y e 27rixy/q 

k=0 Q x:a x =a k 



r-1 -| L(9-fe-l)/rJ 

^^(a'modJV)^ e 2wi(br+k)yq 

k=0 Q 6=0 



The sum taken over all x satisfying a x = a k mod N, or equivalently x = k mod r, is 
replaced by the sum over b, defined by x = br + k, and it is decomposed into two terms. 
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, L(5-fc-l)/rJ 

- E e 

9 fro 



2iri(br+k)y/q 



1 L(«-fc-i)/*-J 

- E < 

? fro 



27rife{rj/} 9 /g _|_ 



|_(,_ fc _l)/ r j 

- E 

? fr 



where {ry} q is congruent to ry mod g (— g/2 < < g/2). The first sum is large 

only when {ry} q is small. When |{ry} g | ~ 0(r), the second sum is 0(1/ q) and the first 
term can be approximated as, 



i/'exp [A] d. + ofi 
r Jo \ r \q 



where the first term is 0(l/r). Neglecting terms 0(1/ q) (<C 0(1 /r)), we obtain the 
probability of observing the state \y), 



r-l 

E e ~ r ~ 

k=0 



iwrr(a k mod N) 



X 



- / exp ( 27ri ^ dw 
r jo \ r 



The second term is the probability of observing the state \y) for a particular value of 
k. If all possible values of k (0, 1, • • • , r — 1) contribute the probability independently, 
the total probability is this term multiplied by r. In our case, the first term in the 
above equation is multiplied instead, due to the fc-dependent phase factors. It can be 
calculated as, 



r-l 

E ' ~ 

k=0 



?-(a k mod N) 



-1 r-l 



= EE 



i^(o* mod N) e -i^a k ' 



e -2N 



k=0 k'=0 

r-l fc-1 



k=l k'=0 



= r + 2 E E cos ^ (( ak mod N ^ ~ ^ mod N )) 



Since —N < (a k mod JV) — (a k ' mod JV) < N, the second term in the above equation is 

k ^ 

positive. The modified visibility X^Z e 1 ^^ mod w ) /r is lower-bounded by 1. 
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